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ABSTRACT 

We investigate the jet propagation and breakout from the stellar progenitor for gamma-ray burst 
(GRB) collapsars by performing 2D relativistic hydrodynamic simulations with analytical analyses. 
We find that the jet opening angle is given by 6j ~ I/SPq, and infer the initial Lorentz factor of the 
jet at the central engine Pq about a few for observations. The jet keeps the Lorentz factor inside the 
star by converging cylindrically via coUimation shocks under the cocoon pressure, and makes a jet- 
breakout acceleration before the free expansion to a hollow-cone structure with a sharp edge. Some 
bursts violating the maximum opening angle 0j,max ^ 1/5 ^ 12° imply a two-component jet. We 
explain both the Amati and Yonetoku spectral relations only with one assumpti on in the off-ce nter 
photosphere model. We also numerically calibrate the pre- breakout model (Bromberg et al.|2011 ) for 
later use. 
Subject headings: gamma rays: bursts — gamma rays: theory — relativity 



1. INTRODUCTION 

Gamma-Ray Burst (GRB) is the brightest object in the 
universe. The observed isotropic energy (the apparent 
energy if it is emitted isotropically) is of the order of 
or even sometimes more than a solar rest mass energy 
MqC^ ~ 2 X 10^^ ergs. Current understanding is that 
the GRB prompt emission is produced by a relativistic 
coUimated jet, whose (half) opening angle is 9j ~ 0.1 
and Lorentz factor is more than P > 100, significantly 
alleviating the energy requirements. 

The opening angle of a GRB jet is an important quan- 
tity not only for the energetics but also for the event rate 
of the GRB. The opening angle also carries information 
of the central engine. It is difficult to get any information 
on the opening angle from an observation of the prompt 
emission, which is beamed into an angle ~ 1/P by a 
relativistic effect. The opening angle of a GRB is mea- 
sured by the light curve of the afterglow that follows the 
prompt GRB. The afterglow light curve exhibits a break 
when the jet decelerates to a Lorentz factor of P ~ 1/^j, 
so-called a jet break, and we can estimate the opening 
angle from the jet break time ( Sari et aL]|1999 ). Recent 
observations suggest that the opening angle ot long GRBs 
are distributed i n several to tens of degrees (6 degrees w 
0.1 radian) (e.g., |Fong eTaLi 2012) . 

What physics determines the opening angle of the GRB 
jet is not known yet. In order to understand the physical 
origin of the jet opening angle, we have to closely look 
into the jet propagation and breakout from the surround- 
ing matter. Because some GRBs are associated with 
supernovae dCalama et al.||1998| llwamoto et al.||1998 



Stanek et al.||20(nf Campana et al. |2006[ Mazzali et al., 
2006p , a long GRB is thought to arise from a death of a 
massive star. A jet is launched deep inside the progeni- 
tor, and should break out from the stellar envelope to be 
observed as a GRB. This is so-called a coUapsar model 



( |Woosley|1993l [MacFadyen fc Woosley|1999[ ) . When the 
jet collides with the stellar envelope, a shocked jet and a 
shocked envelope move sideways from the jet head and 
form a cocoonjj At the expense of the shocked matter, 
the jet head moves outward and finally drills a hole in 
the stellar envelope. This is called the jet breakout |j 

There are many si mulations of the jet propagation for 
the coUapsar model. Aloy et al. (2000) firstly show that 



a relativistic jet can penetrat e the stellar envelope with 



kee ping a good co l hmati on. Zhang et al. (2003 20041 
and Mizuta et al. (2006) study the jet propagation tor 
collapsars with a wide range of jet parameters, s uch as 
the luminosity, th e initial Lorentz factor and so on. |Mor-| 
Sony et al. (2007) discuss the evolution of the jet opening 
angle and indicate that the opening angle of the jet is rel- 
atively smaller than the initial opening angle. However 
the physical origin of the jet opening angle is still unclear. 
What determines the opening angle of the jet? 

One would consider that the opening angle is given by 
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(1) 



where P is the Lore ntz factor at the jet breakout from 



the stellar envelope ( Matzner| |2003 Tom a et al.||2007 ) 
In order to infer P , the cocoon pressure is irnportant as 
it confines the jet ( |Matzner||2003[ |Toma et al.||2007[ lloka 
et al.|20lT|). After passing t hrough the coUimation shock 
(Komissarov & i<'alle 1997), the jet becomes cylindrical 



^ In some articles, a cocoon indicates only the shocked jet which 
has moved sideway. However, the shocked jet is mixed with the 
shocked stellar envelope by the shear interaction and loses its iden- 
tity. Since it is difficult to define the contact discontinuity, we use a 
"cocoon" for both the shocked jet and the shocked stellar envelope 
in a broad sense. 

^ To be precise, a jet breakout is different from a shock breakout 
that occurs when a forward shock reaches the stellar surface. This 
is also different from a jet break seen in the afterglow light curve. 
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since the cocoon pressure is uniform inside the s tar. A 
cyhndrical jet is actually observed in simulations ( Zhan; 






* 



et al.|20 03','2004"Mizuta ct al.|2006[|Morsony et af.|i^OO . 
azzati c t al. 2009, 2013; Mi^ta et al.|201ip and recently 



discussednjy [Broniberg eirar i 2Ull[ ) m an analytic way. 
The Lorentz lactor ot a cyhndrical (stationary) jet is con- 
stant r ~ Fq because of the flux conservation. Therefore 
a naive expectation is that the opening angle is given by 
the inverse of the initial Lorentz factor, 
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(2) 



as shown in the right panel of Fig. nl (conventional pic- 
ture). 

In this paper, we explore the evolution of the opening 
angle for the longGRB jet, and show that the naive ex- 
pectation in Eq. |2| is partly false. We perform a series of 
numerical hydrodynamic simulations for jet propagation 
and breakout from coUapsars. We find that the opening 
angle is a factor ^ 5 smaller than I/Fq, and identify its 
physical origin with the jet-breakout acceleration, that 
is, the jet accelerates at the jet breakout before a free 
expansion as shown in the left panel of Fig. nl (new pic- 
ture). The jet-breakout acceleration boosts the Lorentz 
factor F by a factor of ~ 5 larger than Fq. Thus Eq. (II]) 
is correct but Eq. (l2| is not correct. 

We also examine Fae jet dynamics in an analytical way. 
We first compare the numerical results wi th analytic for - 
mulae before the jet breakout by Bromberg et al. (2011 ). 
Calibrating the model parameters with careful numerical 
calculations, we make the formulae easier to use. Then 
we make analytical analyses for the jet-breakout acceler- 
ation after the jet breakout. 

We further discuss implications of our results. We in- 
fer the initial Lorentz factor of the jet ejected from the 
central engine by using the observed opening angles. We 
also argue a possible origin of the observe d spectral re- 
lations, the Amati and Yonetoku relations (lAmati et al. 
|2002| [Yonetoku et al.| [2004). We derive both the rela- 
tions with only one assumption Lj oc Fq, that the total 
jet luminosity is proportional to the initial Lorentz fac- 
tor, under the photosphere model of the GRB prompt 
emission. These are interesting suggestions for the emis- 
sion mechanism of the prompt emission as well as the jet 
formation mechanism. 

The paper is organized as follows. In Sec. [2] we de- 
scribe the numerical method and the initial conditions 
for the jet parameters. We also introduce the probe par- 
ticles to follow the jet opening angle. In Sec. [Sj we show 
the main results of the hydrodynamic calculations and 
the time evolution of the opening angles. In Sec. H] we 
present the analytic model of the jet dynamics and the 
opening angle of the jet, and compare it with the nu- 
merical results. In Sec. [Sl we apply our results with 
the observations of the opening angle and the spectral 
relations to probe the initial conditions of the jet from 
the central engine and the emission mechanism of the 
prompt emission. Finally, we summarize our results and 
give conclusions of this work in Sec. [6] 

2. NUMERICAL METHOD 
2.1. Numerical scheme 



We have performed two dimensional axi-symmetric rel- 
ativistic hydrodynamic simulations of jet propagation be- 
fore and after the eruption from a progenitor surface in 
order to obtain the final opening angle of the jet dj. An 
updated version of the relativistic hydrodynamic code 
developed by one of the authors (AM) is used for the 
hydrodynamic simulations. The code solves special rela- 
tivistic hydrodynamic equations. We adopt an adiabatic 
equation of state, P = (7 — l)pe, with a constant specific 
heat ratio, 7 — 4/3, where P is the pressure, p is the rest 
mass density, and e is the specific internal energy. 

The hydrodynamic code employs the Godunov type 
fluxes, i.e., an approximate Riemann solver. The versi on 
of Marquina's flux formula ( Donat fc Marquina||1996 l is 
used for numerical fluxes. I'he 2nd or der accuracy in 



space is achieved by MUSCL method (van Leer 
and the 2nd order accuracy in time is achieved 
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total variat ion diminishing (TVD) Runge-Kutta method 
(Shu & Oshcr 1988). Any radiative processes, such as 
emission or absorption of photons are not included in the 
calculations. The details of the code and results of ID 
and 2D test calculations are presented in the Appendix 
of Mizuta et al. (2004, 2006). 

2.2. Grid 

A cylindrical coordinate (z, r) is employed for hydro- 
dynamic simulations (z is the direction of jet propaga- 
tion). In this study, an axi-symmetry is assumed for 
simplicity. 6400 (z) x 500 (r) grid points are used. 
The high resolution grid points are devoted around the 
cylindrical axis so that the interaction with the stellar 
envelope is well captured in good accuracy. The fine 
grid points (10^ cm x 10^ cm) are spaced uniformly at 
lO^cm < z < 4 X 10^° cm and < r < 2 X 10^ cm 
which covers the jet and cocoon at the jet breakout time. 
Logarithmic linear grids are spaced both in z and r co- 
ordinates at 4 X 10^° cm < z < 4.3 x 10^^ cm and 
2 X 10^ cm < r < 1.1 X 10^^ cm. At the outer boundaries, 
the grid sizes are Az = 7x 10^ cm and A?- = 1.6x 10^ cm, 
respectively. 

We perform a resolution study with twice better res- 
olution of the grids [9728 (z) x 500 (r)]. The fine grid 
points (5 X 10^ cm x 5 x 10^ cm) are spaced uniformly at 
lO^cm < z < 4 X 10^° cm and < r < 10^ cm which cov- 
ers the jet and some parts of the cocoon at the jet break- 
out time. Logarithmic linear grids are spaced both in z 
and r coordinate at 4 x 10^" cm < z < 6.2 x 10^^ cm and 
1 X 10^ cm < r < 1.1 X 10^^ cm. At the outer boundaries, 
the grid sizes are Az = 2 x 10^ cm and Ar — 2.2 x 10^ cm, 
respect ively . The results of the resolution study are given 
in Sec. |3Jl 



The resolution with Az„ 



Ar„ 



lO^cm is com- 



parable w ith the highest resolution zone in |Morsony et| 
al. (2007) who adopted the AMR technique, while our 
High resolution area covers much larger area. Our reso- 
lution study uses Azmin = Armin = 5 X lO^cm, which is 
one of the highest resolutions so far. 

The boundary condition is the free boundary condi- 
tion except for the jet injection region and the cylindri- 
cal axis. The reflective boundary condition is imposed 
at the cylindrical axis. 

2.3. Stellar model 
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Fig. 1. — Physical pictures of the jet evolution at the jet breakout. Inside the star, the Lorentz factor of the jet keeps the initial value 
~ To since the jet is almost cylindrical before the jet breakout. After the jet breakout, the opening angle of the jet is ^ I/SFq (left 
panel) rather than ~ I/Fq (right panel), because the jet-breakout acceleration occurs at the jet breakout where the pressure of the cocoon 
decreases outward. 

Assuming the jet formation deep inside the progeni- 
tor, we start the numerical calculation at the distance of 
•Zmin = 10^ cm from the center of the progenitor. From 
the innermost computational boundary we inlet the jet 
into the computational domain. 

At least, four parameters are necessary to characterize 
the initial condition of the jet. We choose the cylindrical 
radius (ro), the luminosity (Lj), the Lorentz factor (Fq), 
and the specific enthalpy {ho = 1-1- eo/c^ -I- Pq/po) of the 
initial jet, where subscription '0' stands for the injection 
parameter. Since we inlet the jet parallel to the jet axis 
with a small radius (rp), we are not necessary to assume 
the initial opening angle of the jet (0o)j which is deter- 
mined by a relativistic beaming effect as da « I/Fq. This 
allows us to reduce the number of the initial parameters. 

(1) We adopt the initial cylindrical radius of the jet 
as ro = 8 X 10^ cm. This value should be sufficiently 




Fig. 2. — Initia.1 rad ial mass density profile (model 16TI in 
[Woosley fc Heger| | |2006[ l). The radial mass density profile is al- 
most power law at lO^cm < iJ < 6 X 10^" cm with the index 
~ —1.5. Then it quickly drops at i? > 6 X 10^ cm. We extend the 
gas to the outside of the progenitor (ij > 4 X 10^'' cm) which is 
assumed to be very dilute and a stellar wind profile with the power 
law index —2. 



smaller than Eq. ( 26 1 which is supposed to be the cylin- 
drical radius of the jet after the collimation shock (i.e., 
the cylindrical radius of the jet balanced with the co- 
coon pressure). If the initial cylindrical radius were 
larger than Eq. (26), the jet dynamics would be differ- 



We adopt one of stellar models in jWoosley fc Heger 
(2006) for the progenitor of the GRB. The model is 



ent from the true one because a thick jet sweeps large 
mass. Even though vq should be small, we still need to 
cover the jet with sufficient grid points. In our simula- 
tions, 8 and 16 grid points cover the jet at the bound- 



named as 16TI. The total mass and radius of the pro- 
genitor at the pre-SN stage are 13.95 solar mass and 
i?* = 4 X 10^° cm, respectively, where i? is a radius in 
the spherical coordinate. The mass density at the inner- 
most boundary (z = 10^ cm) is about 10^ g cm~^. The 
radial mass density distribution of the progenitor is al- 
most power-law with an index ^ —1.5 at 10^ cm < i? < 
6 X 10^ cm and quickly drops at i? > 6 x 10^ cm. The 
mass density at the last grid of the progenitor {R = R^) 
is 1.7 X 10^^ g cm^"^. Outside the progenitor we put a 
low density gas assuming a stellar wind with the power 
law index -2, i.e., p{R) = 1.7 x 10-"(i?/i?*)"2 g ^^-3^ 
See the radial mass density profile in Fig. [2j 

2.4. Jet conditions 



ary for the resolutions. A;; 



= Ar 



= 10' 



and 



Azinin = Armin = 5 X 10 Cm, respectively. Our reso- 
lution is high enough to resolve internal structures, such 
as shocks and vortices, inside the jet and the cocoon. 

(2) We consider a constant luminosity jet {Lj = 5 x 
10^" erg s~^) for simplicity. Thus the explosion energy 
is ^ 10^^ ergs for a few tens seconds duration of the jet 
injection, which is comparable with the values inferred 
by the observations. 

(3) Another two parameters to define the initial jet are 
the specific enthalpy (/iq) and the Lorentz factor of the 
jet (Fq). We fix the product of two parameters, i.e., HoTq 
to be 538 in this paper. The product hoTo gives the maxi- 
mum Lorentz factor achieved by the adiabatic expansion, 
since hT is conserved by the relativistic Bernoulli's prin- 
ciple along a stream line for a steady state. A gas with 
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large enthalpy {h >• 1) expands adiabatically by decreas- 
ing h and increasing F, i.e., the jet accelerates with fixing 
hT. Our assumption {HoTq = 538) satisfies hoTo > 100 
which is required for avoiding the compactness problem 
of the G RB. The recent Fermi bursts suggest a relatively 
large hV ( |Abdo et al!]|2009[ [Ackermann et al.||2010 [ |foka 
2010] ). 

[3] The initial Lorentz factors, Fq = 2.5 , 5, and 10, are 
studied to see the dependence on Fq for the jet dynamics 
and the opening angle after the jet breakout. These are 
the models, G2.5 (Fq = 2.5), G5.0 (Fq = 5), and GIO 
(Fq = 10), respectively. We also perform the hydrody- 
namic simulations with higher resolutions. These are the 
models, G2.5H (Fq = 2.5), G5.0H (Fq = 5), and GIOH 
(Fq = 10), respectively. Since we fix hoTo{= 538), the 
initial enthalpy (ho) is the smallest {ho = 53.8) for the 
case with Fq = 10. Thus all jets are initially thermal- 
dominated plasma {Hq ^ 1). The initial specific internal 
energy (eo/c^) is 80 for the model with Fq = 5. As shown 
later, these initial Lorentz factors are crucial parameters 
for the final opening angles. We set the velocity vector 
of the jet initially parallel to the z-axis. The jet expands 
with an initial opening angle ^ I/Fq as long as the in- 
jection angle is less than ~ I/Fq. 

Table [l] summarizes the jet initial conditions of our 
models. 

2.5. Probe particles 

In order to follow the Lagrange motion of the fiuid ele- 
ments, we introduce probe particles to trace the path of 
the fluid elements. It is necessary to follow the Lagrange 
motion of the fiuid elements to define the opening an- 
gle of the jet, since the jet opening angle depends on 
the time. Every 0.01 s, 32 particles are injected into the 
computational domain with the jet. At the injection re- 
gion at Zmin = 10^ cm, 32 particles are uniformly spaced 
at0<r<ro = 8x 10^ cm. In every hydrodynamic 
time steps (At), the particles move with their local ve- 
locities calculated by the hydrodynamic simulation, i.e, 
^^ncw = iCoid + vAt, where x is the position of a particle 
and V is the local velocity. 

Tracing the Lagrange motion of the fluid elements al- 
lows us to find the location where and how the free ex- 
pansion starts after the jet breakout. Since we define the 
jet open ing angle by the direction of free expansion (see. 
Sec 3.4), the information of the Lagrange motion is very 



important for quantitative analysis, especially to find a 
relation between the opening angle and Lorentz factor. 
The particle path is also useful to identify whether the 
fluid elements have moved into the cocoon or not. 

3. RESULTS 

3.1. Overall evolution 

The jets in all models of Table [T] succeed to drill the 
stellar envelops. Before the jet breakout, the jet hits at 
the reverse shock that is produced by the interaction be- 
tween the jet and the stellar envelope. Then the shocked 
jet moves sideway, forming a high pressure cocoon, as 
shown in the schematic figure |3] The high pressure in 
the cocoon works to confine the jet. The jet is confined 
through a coUimation shock which deflects the velocity 
vector into a parallel direction to the cylindrical axis. 
The high pressure cocoon is also discussed in the context 




Jet head 



Forward (Bow) shock 
Stellar Envelope 

Fig. 3. — A schematic figure of the jet during the forward shock 
is inside the progenitor star. 



of FRII jets dBegelman fc Cioffi|[T989l ) . 

When the forward shock reaches the stellar surface, 
the shocked stellar envelope starts to expand into the 
circumstellar matter, i.e., a shock breakout. Soon the 
jet also start to expand into the circumstellar matter, 
i.e., a jet br eakout as presented in previou s numerical 



simulations (|Aloy et al.||2000| [Zhang et al.|[2b03, 2004 
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the jet breakout, t. 



Morsony et ai.||2d07| [Mizuta fc Aloyj 



'1 



2011; Nagakura et al.||2011[ ). After 



le jet advances in the circumstellar 



matter which is assumed to be very dilute. 

3.2. Before the jet breakout 

The jet is well confined by the high pressure cocoon 
before the jet breakout. Figure [4] shows the contours of 
the mass density, the pressure and the Lorentz factor of 
the jet just before the jet breakout for the model G5.0 
(Fq = 5). We should note that the r axis is elongated 
and the aspect ratio of z and r is not unity in order to 
highlight the fine structures in the jet and the cocoon. 
Although we inlet a jet with a velocity parallel to the 
cylindrical axis, the jet tries to expand with an opening 
angle of 6*0 ~ Fq"^ because the thermal energy is large. 
Around the injection point, the Lorentz factor rapidly 
increases and then drops across a discontinuity. Since 
the pressure and the mass density increase across this 
discontinuity, it is the shock surface. The shock deflects 
the velocity vector almost towards the jet axis, and is 
so-called the collimation shock. The coUimation shock 
is produced by the interaction bet ween the expanding 
jet and the high pressure cocoo n (Komissarov & Falle 
19971 |1998[ [Bromberg fc Levins on 2009; Brombcrg et 



201l| ). To distinguish it from the other coUimating 
oblique shocks in the jet, we sometimes call it the first 
collimation shock. See a schematic figure |3] 
As shown in Fig. |4J after the first collimation shock. 
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TABLE 1 
Jet initial conditions of the models 



model 


To'' 


luminosity (Lj) 
(erg s-l) 


initial jot cylindrical 
radius (ro) (cm) 


hoTo{= Too) 


Az 


mill — ^''"min 

(cm) 


G2.5 
G5.0 
GIO 


2.5 
5.0 
10 


5 X 10^" 


8 X lO'' cm 


538 




10^ '^ 


G2.5H 
G5.0H 
GIOH 


2.5 
5.0 
10 


5 X lO^o 


8 X 10^ cm 


538 




5 X 10" <i 



'^ Initial Lorentz factor 

^ Highest resolution grid size in the computational domain 

'^ at the region of lO^cm < z < 4 X 10^'' cm and < r < 1 X 10^ cm 

'^ at the region of lO^cm < z < 4 X 10^" Vcrmcra and < r < 1 X 10^ cm 



the .jet k eeps an overall cylindrical shape (Bromberg et 
al.|2011 |. The Lorentz factor increases to a tew tens and 
then drops to ^ Fq, i.e., the initial Lorentz factor, after 
the first collimation shock, see Eq 
stand the cylindrical evolution as fo 



(|20|. 
Tows 



We can under- 
The propaga- 



tion velocity of the jet head is comparable to the sound 
velocity of the shocked envelope (obviously) and is much 
smaller than the sound velocity of the shocked jet (with 
high entropy) at the jet head. The shocked jet goes back- 
ward around the jet and provides an constant cocoon 
pressure over the jet. Namely the flow in the cocoon is 
sub-sonic. Although the mean sound velocity of the co- 
coon, which is comparable to the transverse velocity of 
the bow shock, is lower than the jet head velocity, the 
matter is only partially mixed in the cocoon. Thus the 
gas in the cocoon can communicate with each other, re- 
sulting_in a homogeneous pressure profile in the cocoon. 
Figure[5]shows the pressure profile at r = 1.8x lO^cm as a 
function of z. At t = 3 s (before the jet breakout), we can 
see a relatively constant pressure profile over the cocoon. 
Since the jet is confined by the homogeneous pressure 
in the cocoon, the jet keeps a cylindrical structure after 
passing the first collimation shock. The Lorentz factor 
should be constant over the cylindrical jet because of the 
flux conservation. The Lorentz factor keeps ~ Fq, though 
some fluctuations appear in the Lorentz factor due to the 
internal oblique shocks. 

The internal oblique shocks occur in the cylindrical jet 
for several reasons. First reason is that the jet is over- 
deflected to the axis after the flrst collimation shock. The 
jet shrinks the cylindrical radius, resulting in an oblique 
shock. Second reason is similar to that for the first colli- 
mation shock. The shrinking jet becomes over-pressured 
and makes a bounce. This is like the initial expansion at 
the injection and thereby leads to the second collimation 
shock. Such a cycle of the expansion and the collimation 
is repeated several times. Note that the jet has larger 
initial cylindrical radius in the second and later collima- 
tions than in the first o ne. Third reason is t he fiuctuation 
of the cocoon pressure ( Mizuta et al.|[2004 ). The cocoon 
is produced by the the shocked jet and the shocked en- 
velope. The shocked jet at the jet head go sideways and 
also backward under the pressure of the shocked enve- 
lope (Mizuta et al. 2010). Together with the shear mo- 
tion between the shocked jet and the shocked envelope, 
a large vortex and a turbulent structure is formed in the 
cocoon. The turbulence in the cocoon makes a perturba- 
tion at the contact discontinuity between the cocoon and 
the jet, which causes the oblique shocks in the jet. Since 



shocks are tilted to the jet axis, we call these shocks as 
oblique shocks in the jet. 

At the last stage of the propagation in the progenitor, 
we identify a jet-breakout acceleration. Since the den- 
sity distribution of the envelope exponentially drops at 
the stellar surface i?* {z — A x 10^" cm, see Fig. ^, 
the jet head advances rapidly (~ c, the light speed). As 
a result the pressure profile in the cocoon can not keep 
constant around the jet head. The pressure profile in the 
cocoon drops near the jet head. See the one dimensional 
pressure profile in the cocoon in Fig. [5]at i = 4.5 s when 
the jet breakout just occurs. The jet can not keep the 
cylindrical structure near the jet head even after passing 
the oblique shock, since the cocoon pressure drops be- 
fore the jet breakout. The jet expands in a decreasing 
cocoon pressure and makes an acceleration by converting 
the thermal energy into the kinetic energy. The Lorentz 
factor of the jet increases to a few tens even after the 
oblique shock. See the Lorentz factor conto ur n ear the 
jet head in Fig. |4]and the discussions in Sec. |4.2| 

3.3. After the jet breakout 

Soon after the forward shock reaches the stellar surface, 
the jet breakout occurs. The jet and the cocoon (a mix of 
the shocked jet and the shocked stellar envelope) start to 
expand into the circumstellar matter which is assumed 
to be very dilute. The expansion velocity of the cocoon is 
comparable to the sound velocity before the jet breakout 
(~ a few tens percents of the light speed) . So the cocoon 
stays near the stellar surface for ten seconds, providing 
a pressure for the jet confinement. As shown in Fig. [5] 
for t = 5.2 s, the cocoon pressure is decreasing outward^ 
ranging from P ^ const to P ex z~'^ . 

On the other hand, the supersonic jet does not notice 
the cocoon profile until a collimation shock is formed. In 
the star, the jet is repeating a c ycle of th e ex pansion and 
the (over-)collimation (see Sec. 3.2 and|3.4|), and a cer- 



tain expansion starts in the outward-decreasing pressure. 
Note that the expansion is off-center near the stellar sur- 
face, not from the stellar center. The off-center origin 
makes the pressure profile shallower than that shown in 
Fig. [5] with the stellar center origin. Since the pressure 
profile is still shallower than P ex z~^, the jet is col- 
limated but the collimati on is not enough to keep the 
cylindrical radius (see Sec. 3.4). This expansi on le ads to 
the jet-breakout acceleration, as shown in Sec. |4.2| Even 
after the last collimating oblique shock, a certain level of 
confinement continues without forming a shock (or with 
a weak shock), but the jet makes an expansion, leading 
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Lj=5x10°° erg/s, Tg^5, Az„,„=Ar„,„=10^ cm 





Lj=5x10^° erg/s, r(]=5.0, Az„,„=Ar„ 


^=10^cmt=04.5[s] 


8e+09 


'vm^itm 


^ 


4e+09 


m 


Jl 


1 


^ r^ 


^1 


-4e+09 


w- 


tSM 


-8e+09 




JM 



M 



6 
1 5 
4 
3 
2 
1 

-1 
-2 
-3 
-4 



2e+10 
z (cm) 



4e+10 



8e+09 




-8e+09 



Fig. 4. — Mass density in g cni^'' (top), pressure in dyn cm" /c^ 
(middle), and Lorentz factor (bottom) contours of the model Fq = 
5 at t = 4.5 s (model G5.0). Note that the aspect ratio of z and r is 
not unity in order to enhance the fine structures in the jet and the 
cocoon. White dashed line in the Lorentz factor contour indicates 
the initial stellar surface. 



to an additional jet-breakout acceleration (see Sec. 4.2). 
Finally the jet starts to expand freely after the pressure 
profile gets steeper than P ex z~^. As time goes on, the 
steep pressure profile moves inward, while the first coUi- 
mation shock becomes large. Then the jet starts a free 
expansion after crossing the first coUimation shock. 

Figure [6] shows the mass density, the pressure, and the 
Lorentz factor contours of the model G5.0 (Fq = 5) at 
t = 6.4 s. There are many oblique shocks in the jet, not 
only inside the progenitor but also outside the progenitor. 




z(cm)@ r=1.8x10 cm 



Fig. 5. — One dimensional pressure profiles at r = 2 X lO'"' cm 
and t = 3, 4.5, 5.2 s for model G5.0. The profile shows the cocoon 
region. The shocked ambient gas appears as a thin shell at the jet 
head. The pressure profile is almost homogeneous before the jet 
breakout (t = 3 s). When the jet breakout occurs at i = 4.5 s, a 
pressure gradient can be seen around the jet head. After the jet 
breakout t > 5.2 s, the pressure profile in the outer cocoon is about 
p oc z~*. Note that the profile becomes shallower if the origin is 
off-center. 

The oblique shocks outside the star are imprinted before 
the free expansion and expanding in a self-similar way, 
because the internal shocks are not d eveloped in f reely 



accelerating flow without confinement ( loka et al.|2 011). 
The jet advances with a velocity close to the speed of the 
light, whereas the cocoon expands with a sub-relativistic 
speed. 

3.4. Opening angle after the jet breakout 
Since probe particles are introduced in the jet (Sec. 



2.5 1, we can trace the particle path and measure the 
opening angle of the jet. Figure [7] shows the traces of 
particles that are injected at the same time {t — 5 s). 
The path of each particle repeats a cycle of an expansion 
and a coUimation inside the star, and the last collimation 
near the stellar surface is not so strong that the cylin- 
drical radius of the jet becomes large. Just outside the 
progenitor (i? > i?* = 4 x 10^° cm), the particles do not 
yet start a free expansion. At some distances which de- 
pends on the particles, the particle path finally becomes 
straight, i.e., the jet freely expands as discussed in the 
previous section. The position where the particle path 
gets straight depends on when and where the particle is 
injected. The locations to start the free expansion are 
~ 5 X 10"'^° cm for early particles which were near the 
head of the jet at the jet breakout. Then the locations 
to start the free expansion move inward for subsequent 
particles. 

Arrows in Fig. [7] indicate the slopes of the inverse of 
the local Lorentz factor (F"-'^) for the laterally outermost 
particle. The left arrow is F^^ just after the last colli- 
mating oblique shock (z ^ 3.3 x 10^" cm). Although the 
jet behind the shock has higher Lorentz factor than that 
in the star, the arrow is pointing outside the jet opening 
angle. Thus the jet is still confined by the cocoon and 
can not expand freely, drawing a concave particle path 
in Fig. [T] Finally the path becomes straight from the 
base of the right arrow (z ~ 5 x 10^° cm). The direction 
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Fig. 6. — Same as Fig. Hlbut after the jet breakout (at t = 6.4 s, 
model G5.0). 

of the arrow almost coincides with the freely expanding 
direction. This means that the opening angle of the jet 
is determined by the Lorentz factor of the flow when the 
free expansion starts as in Eq. ([T]). 

The extrapolations of the free expansion lines do not 
cross the center of the progenitor, i.e., the explosion is 
off-center. The explosion center moves gradually inward 
as time passes. The off-center position is different even 
for particles injected at the same time. 

Since the particles with hT < 100 lose their potential 
to reach F > 100 by an adiabatic expansion, the only 
particles with hT > 100 contribute to the GRB prompt 
emission. The baryon loading occurs for the particles 



once involved in the cocoon component before the free 
expansion, because the turbulence in the cocoon mixes 
the shocked jet and the shocked stellar envelope. The 
baryon loading also occurs at the contact discontinuity 
between the jet and the cocoon even after the jet break- 
out (mostly via numerical diffusion). A certain level of 
baryon loading is unavoidable through the numerical dif- 
fusi on, a nd we should be careful about it as discussed in 
Sec. 



exhi 



3.7 Most particles which are injected in early time 



)it hr < 100 at large z. 

We measure the jet opening angle {6j) by the angle 
between the jet axis and the free expansion path for the 
laterally outermost particle. Figure[8]shows the time evo- 
lution of the jet opening angle for each model. The solid 
lines show the opening angles measured by the particles 
with hr > 100, while the dashed hues have hT < 100 
that can not produce GRBs. The time axis in Fig. |8] 
corresponds to the time when the particles are injecteoT 
Since the jet breakout time depends on the models, we 
align the shock breakout time at i = in Fig. [8J Thus 
the time t < indicates that the particles move into the 
cocoon before the shock breakout. 

The opening angle of the jet is not constant in early 
times as shown in Fig. |8] There are particles with 
hT < 100 before and a few seconds after the shock break- 
out. Those particles are the components that are en- 
gulfed into the cocoon before the jet breakout. The co- 
coon is largely baryon-loaded because the shocked enve- 
lope mixes with the shocked jet by the shear interaction 
through the contact discontinuity. Since the particles 
expand with the cocoon after the shock breakout, those 
opening angles are relatively large. 

The opening angle drops after the shock breakout in 
Fig. |8] At the same time, the baryon-poor GRB com- 
ponent hr > 100 appears. We note that our results are 
not sensitive to the threshold value of hT > 100 because 
the baryon-loaded flow and the baryon-poor flow are well 
separated, irrespective of a certain un certa inty in the nu- 
merical baryon loading (see also Sec. 3.5). The opening 
angle settles down to a nearly constant value a few sec- 
onds after the shock breakout. Small time variabilities 
of the opening angle are caused by the fluctuations of 
the local Lorentz factor just before the free expansion, 
according to Eq. ([!]). The local Lorentz factor depends 
on where and how the particle crosses the collimation 
shock. Since the size of the first collimation shock grows 
up to a comparable size of the stellar radius, the later 
particles only pass the first collimation shock before the 
free expansion. 

Most importantly, the opening angles in Fig. [8] are 
much smaller than Fq , contrary to our naive expec- 
tations (1/Fn - 0.2 rad for Fq = 5 and I/Fq ~ 0.1 
rad for Fq — 10). Figure^ shows 0j — Fq^ plot for the 
models. We plot the openmg angles at different times 
(from t = 3 s to 7 s) since they fluctuate as shown in 
Fig. [8] The opening angles of the jet are not on the line 
6j ^ Fq ^ but on the line 6j ^ (5Fo)~^. This is our main 
results; the opening angle of the GRB jet from collapsars 
is roughly given by 
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Fig. 7. — Particle trajectories injected with the jet at the same time (i = 5 s, model G5.0). The arrows indicate the slope of the inverse 
of the local Lorentz factor (F"^) at the time when the particle has passed the final coUimating oblique shock before the free expansion 
(left arrow) and at the time when the free expansion starts (right arrow) for the laterally outermost particle. The free expansion direction 
almost coincides with the slope of the inverse of the local Lorentz factor (F~^) when the free expansion starts. The inset displays a zoom 
of the range 1 X lO" cm < 2 < 8 X 10^*^ cm. Note the aspect ratio of z and r is not unity. 



3.5. Angular profile of the jet after the jet breakout 

Figure [lO] shows the angular distribution of isotropic 
luminosity (Liso) using the hydrodynamic quantities of 
the fluid calculations without any constrain for hV, at 
Ra = 1.5 X 10^^ cm at the time when the forward shock 
reaches at zps = 3 x 10^^ cm, for models G2.5 {t = 
14.8 s), G5.0 (t = 12.5), GIO {t = 11.9 s). The angle 
is in spherical coordinate whose origin is the center of 
the progenitor. The angular dist ribu tion of the isotropic 
luminosity profile shown in Fig. [TO] roughly corresponds 
to the particles at t = {zps ~ RaJJc — 5 s after the shock 
breakout in Fig. |10[ The arrows i ndic ate the opening 
angle of the jet measured in Fig. [lO] for each model. 
Each arrow roughly coincides with the rim at which the 
isotropic luminosity starts to drop exponentially. Our 
results are not sensitive to the threshold value of hT > 
100, irrespective of a certain uncertainty in the numerical 
baryo n lo ading. 

Fig. 10 shows a hollow-cone jet structure. The angular 
distribution of the isotropic luminosity is high at the rim 



and drops exponentially at the edge. Then it gradually 
decreases at large angles. The high isotropic luminosity 
rim part is produced by the shock between the expanding 
jet and the high density cocoon before the free expansion. 
The other oblique shocks produced in the star also make 
the structured jet at large z in Fig.[6|and[T0l 

iet are also 



Angular di s tribut i ons of the jet 



Zhang ct aU (20041; Morsony et al.| (|2007[); [Mizuta fc 



TDoy^ (2009tT^ Since they take the radial integration of 



idia 



shown by 



the energy density or the time integration of the energy 
flux at a certain radius, not a snapshot of the isotropic 
luminosity, we can not simply compare their results with 
ours. In fact there are some differences, but it is difficult 
to identify the reasons such as the numerical diffusion, 
the initial jet size, and so on. 

3.6. Fq = 2.5 case 

The behavior of the jet for the case Fq = 2.5 is some- 
what different from other cases (Fq = 5 and 10). The 
jet opening angle at the injection point for the model 
Fq = 2.5 is larger than those for other models, since 
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Fig. 8. — Opening angles as a function of time before and after 
tlie shock breakout for the resolution A^min = Armin = lO^cm 
cases (models G2.5, G5.0, and GIO) are shown. The open ing an- 
gles are measured by using probe particles (see Sec. |3.4[ l. The 
components which satisfies hT > 100 are indicated by sohd lines. 
Those are components to be accelerated F > 100 by the adiabatic 
expansion and can contribute to the prompt emission. The com- 
ponents with /iF < 100 via the baryon loading from the shocked 
stellar envelopes are shown by the dashed lines. 
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Fig. 9. — Opening angle (Bj) — inverse of the initial Lorentz 
factor (Tq ) plot. The opening angles at t = 3, 4, 5, 6 and 
7 s are shown for models G2.5, G5.0 and GIO with the resolution 
^^min = ^''min = 10^ Cm. The Opening angles are also shown at 
t = 3, 4, 5, 6 and 7 s for models G5.0H and GIOH and at t = 5, 6, 
and 7 s for model G2.5H with the resolution Azniin = Armin = 
5 X 10'' cm. We also plot lines of 9j -- 
by naive expectation) and dj = (5Fo)~ 
by numerical simulations). 



Fq (conventional picture 
(our new picture obtained 



the initial opening angle is ^ Fq' . The cylindrical ra- 
dius of the jet becomes larger than those of other models 
(see the analy tic s tudy of the jet dynamics in the pro- 



genitor in Sec. 4.1). As the cylindrical radius of the jet 



increases before the jet breakout, the momentum flux per 
area to push the stellar envelopes decreases, resulting in 
the strong reverse shock. As a result, the forward shock 
and th e re verse shock go away from each other as shown 
in Fig 12 The reverse shock is far from the progenitor 
surface at z ~ 2.4 x 10^° cm at the time when the for- 



FlG. 10. — Angular energy density distribution after the jet 
breakout at a spherical radius {R = 1.5 X 10^^ cm) for models 
G2.5 at t = 14.8 s (green), G5.0 at t = 12.5 (red), and GIO at 
t = 11.9 s (blue). The time is when the forward shock reaches at 
z = 3 X 10^^ cm. The energy density distribution of the jet shows 
a hollow-cone structure. The jet opening angles at i = 5 s in Fig. IS] 
are indicated by the arrows for each model, and coincides with the 
position at which the energy density starts to drop exponentially. 

ward shock reaches the stellar surface, z ^ 4 x 10^° cm. 
Since the cylindrical radius of the jet is very large, a large 
fraction of the shocked jet and the shocked envelope can 
not go sideways and remains in the jet head. The mass 
is collected at the jet head like a snq wplow or a plug 



( [Zhang et al.||2004[ [Mizuta et al.||2006| ). Even after the 
shock breakout, the plug remains on the axis and affects 
the jet advance for a while. As the plug moves away 
from the star, the jet can go around the plug more easily 
and the effect of the plug decreases. A similar struc- 
ture would also appear for th e low luminosity je ts with 



Lj < 5 X 10^" erg s'^ (e.g., Toma et al. 2007), while 



the plug would be reduced for a non-axisymmetric case 
( Zhang et al.|2004 ). We will study the luminosity depen- 
dence in the near fut ure . 

The plug affects the time evolution of the opening angle 
after the jet breakout. The opening angle in the case of 
Fq = 2.5 does not drop quickly, compared with the other 
models. It takes a few seconds to start to drop and gets 
much smaller angle than F^^. As shown in Fig. 12 
there is a baryon-rich plug ahead of the jet when tEe 
jet breakout occurs. The jet is scattered by the plug for 
a few seconds after the jet breakout. As a result the 
opening angle of the jet gets large for a while. When 
the following jets do not interact with the plug which 
moves far away from the progenitor, the opening angle 
gets small as other models. 

3.7. Resolution study 

We have done hydrodynamic sim ulations with higher 
resolution as discussed in Sec. |2.2| The initial jet con- 
tains a small amount of baryon, while the stellar en- 
velope contains abundant baryons, i.e., hT ~ 1. As 
the jet proceeds in the stellar envelope, the contact dis- 
continuity is formed between the shocked jet and the 
shocked stellar envelope. Our numerical calculations use 
the fixed grid points. Artificial baryon loading could hap- 
pen when the discontinuity crosses the grids via the nu- 
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Fig. 11. — One dimensional profile of the Bernoulli's constant hr 
along the jet axis for Fq = 5 at t = 6.4 s when the forward shock 
reaches z ~ 10^^ cm. Both resolution cases, i.e., Azn^in = Arniin = 
lO^cm (model G5.0) and Azmin = ^''min = 5 X lO^cm (model 
G5.0H) are shown. hF is conserved up to the jet head except some 
fluctuations at the internal shocks. On the other hand, hF drops 
near the jet head for a simulation with similar jet parame ters but 
a low-resolu tion spherical coordinate, A9 = 0.25 degrees, | |Mizuta| 
|et al.|2011| l due to the baryon loading via the interaction with the 
stellar envelope. 



merical diffusion. Once the gas is polluted by baryons 
{hT ~ 1), it is not accelerated to a large Lorentz fac- 
tor by the adiabatic expansion. Figure [TT] shows one 
dimensional profile of hT along the jet axis for Fq = 5 at 
t — 6.4 s when the forward shock reaches z ~ 10^"*^ cm 
with both resolution cases, i.e., Az,nin — Ar„iin — lO^cm 
and Aznun = Armin = 5 X lO^cm. In both cases, the 
Bernoulli's constant hT is conserved all the way to the jet 
head except for some fluctuations at the internal shocks. 
If we apply poor resolution, the Bernoulli's cons tant hT 
is not well c onserved. For example, the case of |Mizuta| 



et al. (2011) is also shown. hT is conserved up to a half 
ot the jet but n ot near the jet head d ue to the numerical 



baryon loading. Mizuta et al. (20111 adopted similar jet 
parameters but a spherical coordinate for hydrodynamic 
simulations with A9 — 0.25 degrees around the jet axis. 
This is one of the reasons to use much higher resolution 
grid points for detailed quantitative discussions on the 
openingangle of the jet. 

Fig. [9| shows the opening angles for different resolu- 
tions. The results are similar between different resolu- 
tions. Therefore the relation 9j ^ (SFp)"^ in Eq. (|3| 
seems robust. 

Although the differences are small between the models 
G2.5 and G2.5II (Fq = 2.5), they depend on the resolu- 
tion most sensitively than the other models. Figure 13 
shows the mass density, the pressure, and the Lorentz 
factor contours of the model Fq = 2.5 with the high res- 
olution at i = 7.5 s when the shock breakout occurs. 
The shock breakout times (6.6 sec and 7.5 sec) are dif- 
ferent by ^ 1 sec even if the jet parameters are the same 

The high resolution case takes 



in Fig. 12 and and 13 



about 1.2 times longer time than the lower resolution 
case. One can notice that the recollimation shocks are 
different. The low resolution one converges on the axis 
while the high resolution one makes a Mach disk. Prob- 
ably the difference is caused by the nonlinear evolution 
of the jet and the cocoon dynamics. A small difference 
in the oblique shocks near the head of the jet changes 
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Fig. 12.— Same as Fig. [I] but Fq = 2.5 (model G2.5), just 
before the shock breakout (i = 6.6 s). 

the cross-sectional radius of the reverse shock and hence 
the jet head speed. This affects the difference in the ap- 
pearance of vortice s and the turbulence in the cocoon 
( Mizuta et al.||201^ ) . The turbulence itself is also a non- 
linear process! The numerical diffusion of the baryon 
loading would also alter the dynamics. The size of the 
head plug is also larger for the high resolution than that 
for the low resolution when the jet breakout occurs. 

3.8. Initial jet size 

At the injection point Zmin = 10^ cm, we set the initial 
cylindrical radius of the jet to be ro = 8 x 10^ cm. The 
initial cylin drical radius should be smaller than the value 
in Eq. ( 26 ) that will be the cylindrical radius of the jet af- 
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Fig. 13. — Same as Fig. [4]but Fq = 2.5 with a high resolution 
(model G2.5H), just before the shock breakout (i = 7.5 s). 

ter the coUimation shock. Otherwise, the jet evolution is 
different. The jet does not show an initial expansion but 
keeps the incorrect cylindrical radius. This gives wrong 
dynamics for such as the breakout time, the cocoon ener- 
getics, and the opening angle. That is why we start the 
numerical simulation with so small initial cylindrical ra- 
dius. At least several or ten grid points are necessary for 
covering the initial cylindrical radius because we have to 
resolve the internal shocks. This requests us to use fine 
grid points as used in this study. Our simulations use 
one of the finest resolutions so far. 

3.9. Comparison with previous studies 



We have done one of the most highest resolution cal- 
culations for hydrodynamic simulations of the jet prop- 
agation. High resolutions reduce the numerical baryon 
loading which affects the jet dynamics. Since the numer- 
ical baryon loading is the most dangerous at the bound- 
ary with a high contrast of the enthalpy h, the highest 
resolution grid points are devoted to the region which 
covers the jet and some part of the cocoon. The high- 
est reso lution grid size is comparable or better tha n that 



used in Morsony et al. (2007); Lazzati et al. (2013). The 



hi ghest resolution region is much larger than that used 



Morsony et al. 



( |2007[ ). 

Since we have done quite high resolution calculations 
for hydrodynamic simulations of the jet dynamics, the 
computational domain is restricted to about ten times 
larger than the progenitor size . This is smaller than 
that used in Mizuta et al. (2011); Nagakura et al. (2011 ); 
Suzuki fc Shigeyama_ (2013) . Our discussion is restricted 
only tor several seconds after the jet breakout. The com- 
putational box si ze is comparable with that used in Mor- 
sony et al. ( 2007 ) who discussed the time evolution of the 
opening angle or the jet. The requirement that the initial 
cylindrical radius of the jet should be sufficiently small is 
also one of th e reasons to use high resolution grid points 
(see Sec. 3.8). We have paid a special care with this by 
taking a sufficiently small cylindrical radius of the initial 
jet. The dyn amics of our results m ay be different from 

((20071 



those done by 



Morsony et al. 



who take the initial 
cm for 6*0 = 10 degrees 



cylindrical radius to ~ 2 x 10 

case. This is about three times larger than our initial 

cyfindrical radius of the jet (ro = 8 x 10^ cm) . 

Our numerical results show that the opening angle of 
the jet after jet breakout is narrower than that we expect 
{0j ^ 1/ro). This trend is consiste nt with previous nu - 
merical simulations, for example, iMorsony et al. (|2007| 



Lazzati et al. (2009); Mizuta et aEJ l 2011)7 Lazzati et al.^ 
( |2013p who have pointed out that tne opening angle after 
the jet breakout is smaller than the initial opening angle 
(^o)- In this paper we introduce much more sophisticated 
analysis with probe particles which allow us to follow the 
Lagrange motion of each element. We find out the corre- 
lation between the initial Lorentz factor and the opening 
angle of the jet, i.e., 9j ^ I/SFq. We also find where 
and how the gas starts the free expansion. The position 
to start the free expansion is off-center and also differ- 
ent for each particle around ^ 4 x lO^*' cm. We identify 
that the jet breakout acceleration occurs and the local 
Lorentz factor of the flow just before the free expansion 
determines the opening angle of the jet. Since we mea- 
sure the opening angle by the free expansion direction 
respect to the jet axis, the angle is different from that 
measured by Morsony ct al. (2007) who defined the jet 
opening angle by measuring how the GRB jet component 
(Finf = hT > Tcr) spreads out at i? = 1.2 x 10^^ cm in the 
spherical coordinate (the progenitor center is the origin 
of the coordinate) . The different measurement methods 
result in about 30 percents difference in the opening an- 
gle of the jet. 

4. ANALYTIC MODEL 

Let us consider the analytic models for the jet propa- 
gation and the dynamics of the opening angle. The evo- 
lution is generally divided into two phases, i.e., before 
and after the jet breakout. 
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In Sec. |4.1[ we discuss the pre-breakout phase. 



Bromberg et al. ( |2011 ) investigated the pre-breakout evo- 
lution in detail. Betore the breakout, the jet is inside the 
progenitor star and cohimated by the cocoon pressure. 
Since the cocoon pressure is constant (subsonic inside), 
the jet becomes cyUndrical through the first coUimation 
shock. The cyhndrical (and stationary) flow has a con- 
stant Lorentz factor because of the flux conservation. If 
the Lorentz factor were constant even at the jet break- 
out, the opening angle of the jet would be an inverse of 
the initial Lorentz factor (or the initial opening angle; 
see below) as shown in Fig . Ill but this is not the case as 
shown in the next Sec. 14.21 

We present analytical formulae of the jet propagation 
for easy comparison with our numerical simulations, and 
make a calibration of the model parameter taking great 
care of the baryon contamination and the initial jet size. 
We compare the evolution of three physical quantities, 
the jet head position, the jet cylindrical radius, and the 
coUimation position, with numerical calculations. The 
model has one free parameter rj in Eq. ([8]), which is ob- 
tained by a fit ting to the numerical results. 

In Sec. |4.2[ we argue the post-breakout phase. Af- 
ter the breakout, the cocoon pressure is not constant 
but decreasing outward (see Fig. [5]). This leads to less 
coUimation, wider jet, and hence larger Lorentz factor. 
Therefore, the jet-breakout acceleration occurs almost 
inevitably if the external medium has a finite size. We 
estimate that the jet-breakout acceleration boosts the 
Lorentz factor of the jet by a factor of several ^ 5. This 
is the reason why the naive picture in Fig. IT] is not cor- 
rect but the jet opening angle becomes ^ l/oT as shown 
numerically in Sec. |3.4| and Fig. |8] 

4.1. Jet evolution in constant external pressure 

Fir st we consider the let evo lution inside the progenitor 
star. Bromberg et al. (2011) give detail analyses. Here 
we make the analytical formulae easier to use in the cal- 
ibration of the model parameters. The jet dynamics is 
controlled by three processes, (1) jet head, (2) cocoon, 
and (3) coUimation. After combining these dynamic s, w e 
can give the evolution of the jet head position in Eq. ( 25 ) , 
the jet cylindric al ra dius in Eq. (26), and the coUimation 
position in Eq. (27 1. There is one free parameter 77 in 
Eq. ([8]) to be fixed by numerical simulations. 

(1) Jet head dynamics: After propagating inside the 
star, the jet collides with the stellar envelope. A re- 
verse shock decelerates the jet, and a forward shock runs 
into the stellar envelope. The shocked region is called 
the jet head. The jet head dynamics is determined by 
the ram pressure balance between the shocked jet and 
the shocked envelope, both of which are given by the 
pre-shock quantities through the shock jump conditions 
( Marti et aLl|1997[ [Meszaros fc Waxman||2001[ [Matzner 



Wm^ , 



h,PjC^r%(3l + P, = KpaC^Tlpl + Pa, 



(4) 



where Tjh — TjTh{^ — PjPh) is the relative Lorentz 
factor between the jet and the jet head, and /3j/i = 
{Pj — /3/i)/(l — I3j/3h) is the corresponding relative ve- 
locity. We can neglect the internal pressure of the jet 
Pj for the strong reverse shock and the pressure Pa for 
the cold ambient matter. Then the jet head velocity is 



obtained as 



where 



/?h 



/3. 



l+i-l/2' 



L: 



ha Pa ^jPaC^ 



(5) 



(6) 



is the ratio of the energy density between the jet and the 
ambient medium. In the last equality, we assume the cold 
ambient medium ha = 1, and use the jet cross-section 
Y}j = nr'^ and the jet luminosity Lj. For typical parame- 
ters of GRBs, we have L ^ 1, i.e., a non-rclativistic head 
velocity, 

Ph ^ i'/'. (7) 

Hereafter we consider the non-relativistic case. 

(2) Cocoon: The shocked jet and the shocked envelope 
try to expand and go sideways into a cocoon component. 
The cocoon pressure is determi ned by the inject ed energy 
divided by the cocoon volume ( [Begelman fc Cioffi_1989| ), 



E ^T] jLj{l-Ph)dt 
3Vc 3iJPhcdt)7T{J/3ecdty 



(8) 



where rj is 0(1) parameter to correct the approximation 
of the cylindrical cocoon shape. We will use rj to ab- 
sorb other approximations, such as that we represent the 
transverse velocity by a single value, and determine 77 by 
comparing the analytic formulae with the numerical sim- 
ulations. Note 1 — /?ft « 1 for the non-relativistic head 
velocity. The transverse velocity of the cocoon is deter- 
mined by the balance between the cocoon pressure and 
the ram pressure of the ambient medium. 



/3c = 



'A. 

PaC^ 



where 



_ / PadV _ 

Pa(Zh) = 77 = UPa(Zh), 



(9) 



(10) 



is the mean density of the medium. Defining 

f /3^dt = ^hf3ht[^^), I pcdt = ^,P,t, (11) 
we can eliminate /3c from Eqs. (l8| and ^ to obtain 

Pr = L-^' 



/ T \ 1/2 / c \ 1/2 



Ct^ 



avr^Cc' 



(12) 



If the density profile follows a power-law pa on z " , the 
coefficients, ^a, £,h and ^c, become constant. 



6 — a 6 



(13) 



Here we obtain ^^ in Eq. ( 10 ) assuming a spherical co- 
coon with a radius Zh for simplicity. Although this is 
not a good approximation, we adjust 77 in Eq. (p| to fit 
the analytic mod els with numer ical results. For ^h and 
^c, we use Eqs. ([22|), ^ and ([23]), which yield /3h oc 

t-^/^Pht)"^^ oc t^ and Pc ex i-2/5(/3?,i)"/5 ex t^ . 
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.'?i=0 



?l=1.{ 



determines the cross-section of the jet after the coUinia- 
tion shock, 



Sj(z > z) 



L, 



70 2n 



4cR, 



(18) 



Combining with a general 
(radiation-dominated jet), 



relation for a hot jet 



L, ~ 4P,r2E^c, 



(19) 



we obtain the Lo rentz factor of the jet after the first 
collimation shock (Bromberg et al.||2011), 



Fig. 14. — Structure of the jet and the coUimation shock. The 
coUimation shock appears by the interaction between the expand- 
ing jet and the high pressure cocoon. The coUimation shock con- 
verges at z. Sohd lines show the colhmation shock and the jet 
structure after the coUimation shock for a decreasing pressure case 
(A = 1.8), whereas green colored dashed lines show the coUimation 
shock and the cylindrical jet structure after the collimation shock 
for a constant pressure case (A = 0). 

(3) Collimation shock: If the ambient density is suffi- 
ciently high like in the stellar envelope, the cocoon pres- 
sure becomes high enough to coUimate the jet. The ini- 
tially expanding jet hits the first collimation shock and 
converges its trajectory. For a constant cocoon pressure, 
the conical jet becomes cylindrical after the collimation. 
This process determines the cross-section of the jet and 
thereby the jet Lorentz factor after the shock. 

The geometry of the collimation shock is determined 
by the pressure balance between the jet and the cocoon 
( [Komissarov fc Falle||1997[ [Bromberg et al.||2011[ ), 

(14) 



hoPoc'T'^/3^ sin^ * + Po = Pc 

where the first term is the ra m pr essure of the jet normal 
to the shock surface (see Fig. 14 1, the subscript stands 
for the unshocked jet. Since the jet internal pressure 
decreases Pq cx z~^ as the size grows, we neglect the 
term Pq. For a small incident angle, we have a relation 



sin * = — - —^ 

z dz 



= z— — 



dz 



m 



(15) 



to a first order (see Fig. 14), where r^ is the cyl indrical 



radius of the shock position. Then equation ( 14 1 gives a 
first-order ordinary differential equation. Notmg /3o ~ 1 
and Lj ~ hopoc^rQ{Trz^9Q), we can integrate the geome- 
try of the collimation shock as 

Ts = Ooz{l + Az^ 



%Az' 



(16) 



where A is given by Eq. (17 1. We note that we assume 
a constant external pressure Pc =const. For a decreas- 
ing external pressure, the shock geometry is different as 
shown in the next Sec. [4.21 

According to Eq. ([16[), the collimation shock expands 
to a maximum size at -p^L_. — 0, and converges to 
rs{z = z) = where the maximally-expanding position 
■Zmax and the converging position z are given by 



2z„ 



A i + z. 



A- 



L,/3, 



jPO 



TTcPr 



(17) 



In the second last equality we assume that the collima- 
tion shock is initially small z^, ^ A"^. As we can see 
from Fig. [14[ the maximum size of the collimation shock 



Ti 



(20) 



After the collimation shock, the jet is usually over- 
defiected to the axis, resulting in an oblique shock inside 
the jet. A high pressure region appears after the converg- 
ing point of the first collimation shock and then expands 
again. The jet repeats a cycle of a bounce, expansion 
and collimation, like a spring under the pressure of the 
cocoon. Because the supersonic jet is not synchronized 
with the cocoon, the oblique shocks arise in the jet. In 
each collimation, the jet tries to expand with the same 
angle ~ I/Fq and hence the Lorentz factor after a colli- 
mating oblique shock keeps the same value. The Lorentz 
factor after the A^-th collimating oblique shock is 



r 



N 



T 



(w-i) 



Ti 



1 



To, 



(21) 



under a constant pressure inside a star. The last equal- 
ity is satisfied in our simulation because we inject a jet 
parallel to the z-axis and the jet tries to expand with an 
angle Oq ^ l/Fg about an inverse of the initial Lorentz 
factor. 

If the Lorentz factor were constant F ~ Fq even at the 
jet breakout, the opening angle of the jet would be the 
inverse of the Lorentz factor inside a star (or the initial 
opening angle) 6j ^ T^^ ^ 9q. Howeve r th is is not the 
case as discussed in the next subsection 14.21 

For comparison with the numerical results, we express 
physical quantities by basic parameters, i.e., the jet lu- 
minosity Lj, the ambient density pa, the initial opening 
angle Oq, and time t. First, we rewrite the head velocity 
and the cocoon pressure as 




L, 



t^PaOt 






1/5 



,3 r2fl2\ ^/^ 



1 vL 



2/5 



(22) 



(23) 



respectively, with equations (|6|, Q, (12 1 and (18). Then 
we can der ive analyt ic fo rmulae iSr the jet head position 
with Eqs . (11) and (22), the jet cylindrical radius with 
Eqs. ( 18 ) ana ( 12 ), ana the converging point of the col- 
limation shock with Eqs. (17) and (12) as 



Zh = ihPhct 



t^L.V fl&vUt^'" 



PaOt 



3^ ec 



(24) 
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(25) 



2\ 1/5 



(26) 






1/10 



!4.7x 10^ cm 



Pa 
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2/5 



2\ 1/5 



i. 



Is/ \ 10^^ ergs" 

-3/10 / g, X -1/5 

01 



3/10 



(27) 



respect ively , where we set a = 2 [^^ — 3, 61 = Cc = 1 
in Eq. (13l] and 77 = 0.01 for the numerical values. For 
a general density profile rather than a power-law form, 
we solve the first-line equation ( 24 ) for the head position 
Zh, where the quantities pa, £,a7Xh a-nd Cc in the right- 
hand-side are also functions oizh- Here, for simplicity, 
we set £,a, £.h and fc with Eq. (13) by using the density 
slope a = —dlnpa/dhiz at z = zu- 

As show n in Sec. [31 we can fit three analytic formulae in 
Eqs. (25), (26) anapTJ) with the numerical calculations 
by adjusting one parameter 77 ^ 0.01 in Eq. ([s]). Figure 



|15| shows the comparison with analytic formulae and the 
results of the hydrodynamic simulations, i.e., the head 
position, the converging position, and the jet radius, for 



the model G5.0 (Fq = 5 and Azj, 



Arn 



lO'^cm). 



The numerical results provide good agreement with the 
analytic ones. The difference in the early evolution of the 
jet is caused by the insufficient smallness of the initial jet 
size. Since the power law index (a) in the mass density 
profile exceeds a > 3 at a certain radius and ^a becomes 
infinity or negative in Eq. (13), we cannot apply the 
analytic formula after t = l.(J5 s. Figure [TS] is the confir- 
ma tion th at the analytic physical picture by Bromberg et 



al. (2011) is correclrjby the numerical calculations that 
take great care of tne baryon loading and the initial jet 
size. 

4.2. Jet evolution in decreasing external pressure 

We numerically find that the opening angle of the jet 
is not ~ I/Fq but ~ l/5Fo where Fq is the Lorentz factor 
of the jet inside a star. This is because the cocoon ex- 
pands to the outer space and hence the cocoon pressure 
decreases outward after the jet breakout from the stellar 
surface. In the decreasing pressure, the jet becomes less 
coUimated and hence more accelerated than that in the 
star. Since the jet-breakout acceleration increases the 
Lorentz factor F before a free expansion, the jet opening 
angle determined by ^ 1/F is narrower than the naive 
estimate. 

^ Althoug h ri ~ 0.01 is smaller than 0(1), a factor ~ 2 difference 
in 2ft in Eq. (|24|) may cause a factor 2^ ^ 32 difference in tj. 
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ical formulae for the jet propagation inside the star (see Sec. |4.1[ l. 
We show the forward shock position of the jet head 2^, the converg- 
ing position of the collimation shock z, and the maximum cylin- 
drical radius of the jet rs^max — ^oZmax- Theoretical model lines 
are also plotted. Since the slope of the radial mass density profile 
starts to rapidly drop at _R ~ 6 X 10^ cm, there is no analytic 
solution outside this point. 

The jet-breakout acceleration is broadly divided into 
two types, i.e., with and without a shock. The both types 
contribute equally to the Lorentz factor enhancement. 

First we consider the jet-breakout acceleration by eval- 
uating the Lorentz factor after the coUimating oblique 
shock. The cocoon pressure that is decreasing outward 
is expressed by 



P,{z) = n ( - 



(28) 



Under this pressure, we can solve the geomet ry o f the 
collimation sh ock , as in the previous section |4.1| with 
Eqs. pll) and ((iSl), 



iz 1 



A,z, 






(29) 



where z* is the initial position of the collimation shock 
and 



A^ — 






(30) 



We show the shock geometry in Fig. [Ml Note that the 
shock geometry is not parabola but slee k in contrast to 
the constant pressure case in Eq. (16). The collima- 
tion shock expands to a maximum cylindrical radius at 
-2^|z=Zniax — 0, where the maximally-expanding position 
is given by 



2-A 



1 



A 



(31) 



Here A < 2 is necessary for the shock to converge. At 
the posi tion Zmax, the general relation for a hot jet in 
Eq. ( 19 ) is given by 



P* 



-A' 



r? {< 



2 2 \ 
in ax/ ^" 



(32) 
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This yields the Lorentz factor after the collimation 
shock, 



Ti 



X A, 



(33) 



which is larger than that for a constant pressure case in 
Eq. pOl) by 



A 



4- A 



2A + 4A. 



4- A 
4-2A 



5, if A - 1. 



(34) 

for small A*z* . To be precise, 6q is not the initial opening 
angle of the jet here, but the opening angle of the jet 
expanding to the last coUimating oblique shock, which is 
the inverse of the Lorentz factor insid e a jet ^ r,^ and 



thereby turns out to be ^o from Eqs. (20) and (21). 

The numerical calculations show that the totaTacccl- 
eration factor is .4 ~ 5 (see Fig. IIs]) , about half of which 
is achieved at the collimation shock and the other half 
of which is obtained later. A factor ^ 2.5 can be ex- 
plained by A ~ 1.5. Although it is difficult to estimate 
the exact acceleration factor A analytically, we can see 
that A is mainly determined by the slope A of the exter- 
nal pressure profile. In our case, the external pressure 
profile is shaped by the cocoon expansion to the outer 
space, which does not depend on the jet properties so 
much. Therefore it is natural that similar acceleration 
factors are obtained for different initial conditions in our 
numerical calculations. 

Next we consider a jet-breakout acceleration without 
shocks (i.e., an adiabatic jet) as the other extreme. The 
jet expands in a decreasing pressure, lowering its temper- 
ature as radiation T' ex V'~^^"^ where V is the comoving 
volume of the jet. If the pressure balances between the 
jet and the cocoon T'* (x Pc, the jet does not perform 
work and hence the energy is conserved TT'^V' ^ const. 
Th erefore the jet ac celerates as the size grows according 
to dloka et al.|l2011|) 



r ex T'~i 



p-1/4 



,\/i 



(35) 



In the numerical calculations, we identify such an adia- 
batic evolution after the collimation shock, in particular 
at the periphery of the jet. For 2 < A < 4, the cocoon 



pressure does not causally affect the jet interior (loka et 



al.||2011[ ) but still affect the periphery of the jet"^ Note 
that the off-center origin makes the pressure profile shal- 
lower than that shown in Fig. [5] with the stellar center 
origin. 

In order to estimate the acceleration factor of the 
Lorentz factor, we need to know how much the jet ex- 
pands before entering a free expansion phase. In this 
regard, we should note that the jet expands from the 
breakout position, that is, the fireball of the jet is off- 
centered by the cocoon confinement. This means that 
the fireball size should not be measured from the center 
of the star. Instead, the effective center of the fireball 



locates at a distance ~ ToVj 



j/00 inward from the 



stellar surface (breakout point) because the jet tries to 
expand with an opening angle of I/Fq and the expanding 
surface has an initial cylindrical radius of rj in Eq. ( 26 ) . 
So the initial fireball size is 



ro ~ Tor J ~ rj/0o. 



(36) 



0) 

E 




Fig. 16. — Distribution of the initial Lorentz factor of the jet 
inferred by tlie observations and tiie relation Fq = 1/50 , in Eq. ||3ll. 
Observational data {9j , N) are taken from Fig. 6 in |Fong et"al.| 
| |2012| . Since a Lorentz factor should be greater than unity, the 
region with Fq < 1 is shaded as an unphysical zone. 

If the fireball expands to a size of the stellar radius 
i? ~ 4 X 10^° cm, the expansion factor is about ~ 10 
times for typical parameters in Eq. ( 26 ) and hence the 
Lorentz factor grows by a factor ^ 3 tor A ~ 2 according 
to Eq. (35). Thus the adiabatic expansion can explain a 
part of the jet-breakout acceleration observed in the nu- 
merical calculations. Here a parameter dependence of the 
expansion factor is weak since the initial fireball size is 

ro ex L Pa 0Q at the breakout time (when R ~ Zh) 



with Eqs. (25| and (26). Note that the jet-breakout ac- 



celeration looks very rapid at the first glance if we do 
not notice the off-center effect (i.e., the radius measured 
from the stellar center is only doubled). 

5. DISCUSSIONS 

5.1. Lorentz factor of the jet in a star 

The opening angle of the GRB jet is usually mea- 
sured by observing a jet break in the afterglow lightcurve 
( [Racusin et "aL||2009[ |20TT| |Fong et al.l|2012[ ). Our nu- 

merical and analytic calculations show that a jet opening 
angle is related with the Lorentz factor inside a star as 



1 



(37) 



By applying this formula, we can infer the Lorentz fac- 
tor inside a star (or the initial opening angle) from the 
observed opening angle of the GRB jet. Figure [T6| shows 
the estimated Lorentz factor inside a star. The result 
suggests that the jet is mildly relativistic in a star (or 
the initial opening angle is 0(0.3-0.5) radian). 

The slow jet inside a star is a preferable condition for 
the survival of nuclei in the jet, which may explain the 
observed u ltra-high energy cosmic-ray nuclei (|Mu rase et 
al. 2008; Wang et al.|2008||Metzger et al.|201in ni 
|et al.||2012| 



loriucni 



In 



Equation ( 37 ) also implies that the maximum opening 
angle is obtamed by setting Fq ^ 1 as 



1/5 = 0.2- 12°, 



(38) 



if the jet is radiation-dominated at the breakout. How- 
ever several events violate this maximum as shown in 
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Fig. [I6| pacusin et al.|[2009l [20Tl| |Fong et al.|[20l2| . A 
possiEIe solution to this problem is to consider a baryon- 
rich slow sheath surrounding a central jet. A baryon-rich 
flow cannot accelerate to a large Lorentz factor, and if 
the Lorentz factor of the baryon-rich sheath is less than 
'^ 5, the opening angle of the sheath can be larger than 
~ 0.2. Note that the central jet should have F > 100 to 
avoid the compactness problem. 

5.2. Amati and Yonetoku relations 

Based on our numerical and analytical analyses, we can 
evaluate the initial condition of a jet that is just expand- 
ing freely, and infer the observational characteri stics of 
the je t based on the photosphere model (e.g., Pe'er et al. 
|2007 loka et al. 2007|. First, since we are now able to 
estimate the opening angle of the jet in Eq. (371, we can 
assess the isotropic luminosity of the jet. 



-^iso — 



L, 



ex i,F| 



■j^ 0- 



(39) 



We can also obtain the observed temperature of a jet if 
it is radiation-dominated. 



E. 



peak 



Tobs (X FoT' (X ToP}'\ 



(40) 



which may be identified with the spectral peak energy 
Epcak of a GRB in the photosphere model. In the actual 
observations of GRBs, there is a relation between the 
isotropic luminosity and the observed temperature. 



1 /? 

Epcak oc L-J^ (Yonetoku relation), 



(41) 



so-called the Yonetoku relation ( [Yonetoku et al.|[2004[ ). 

Let us show that the Yonetoku relation is reproduced 
if the total jet luminosity is nearly proportional to the 
Lorentz factor inside a star 



Lj oc Fo, 



(42) 

that is, the mass outflow rate is M = Lj/Tq ^ const. 
First, the above equation give s Liso oc Fq with Eq. (39). 
Next, substituting Eqs. p9| and (|40| into Eq. (T9| 
(where we should read Fi as Fq) yields 



p ^ 7-1/4^-1/4 



(43) 



Here the jet cross-section follows Ej oc L-' 9q' at the 
time 



t oc L, 



-l/3/)4/3 



(44) 



of the jet breakout, Zh ^ R* (when the jet head reaches 
the stellar surface) with Eqs. (25) and ([26|F] Then, not- 
ing Oq ^ Tq , we have 

E, 



rl/4p7/12 .4/9 
OC iv:„„ in OC iv- 



(45) 



which is close to the Yonetoku relation. 

In additio n, we may be also a ble to reproduce the Am 
ati relation ( Amati et al.||2002 ). 



1 In 

^pcak ^ E-Jo (Amati relation) 



(46) 



4 We assume that the structure of the progenitor star is similar 
and neglect the dependence on it. 



We can think that the GRB duration is roughly given by 
the sound-crossing time across the cocoon of the stellar 
size, 



Tc 



90 






t^/^L] 



-l/5/)-l/5 



OC L, 



-l/3/)l/3 



« L7.rm 



with Eq. ^ and (23), because the cocoon pressure de- 
creases after this time, leading to less confinement, larger 
opening angle, and smaller isotropic luminosity of the jet. 
Note that the weak correlation between Tgp and Ljgo is 
actually observed. Then we can estimate the isotropic 
energy 



E, 



peak 



oc Lfi' oc (Li^oTgo)"/' ^ E. 



4/7 



(48) 



which is also similar to the Amati relation. This is killing 
two birds with one stone, that is, we explain two relations 
(Amati an d Yonetoku relations) by only one assumption 
in Eq. ^. 



6. SUMMARY 

In this paper we have explored the dynamics of the 
GRB jets from coUapsars by performing 2D relativis- 
tic hydrodynamic simulations as well as the analytical 
analyses. We have followed the jet propagation all the 
way from the progenitor star through the jet breakout to 
the free expansion, with implementing probe particles to 
trace the Lagrange motion of the fluid elements. This en- 
able us to connect the final jet appearance to the initial 
jet condition near the central engine. 

We have found that the jet opening angle after the 
jet breakout is about Oj ~ l/5Fo in Eq. ([sl and Fig. 
where Fq is the initial Lorentz factor of the jet injecte' 
into the progenitor star. This value is smaller than the 
naive expectation 9j ~ I/Fq in Fig. [11 where we thought 
that the opening angle is determined by the inverse of 
the Lorentz factor and the Lorentz factor keeps the ini- 
tial value Fq for a cylindrical jet. Actually this is partly 
correct. The jet becomes cylindrical under the nearly ho- 
mogeneous pressure of the cocoon after crossing the first 
coUimation shock. The Lorentz factor after the coUima- 
tion shock is ~ Fq and largely stays constant before the 
jet breakout according to our simulations. However we 
have identified the "jet-breakout acceleration" (just be- 
fore and) after the jet breakout. This occurs because the 
pressure profile of the cocoon can not keep constant but 
decreases outward as the cocoon expands to the outer 
space. The cocoon still confines the jet to some extent 
near the stellar radius, while the jet expands its cylindri- 
cal radius with increasing its Lorentz factor by a factor 
'^ 5 before a free expansion. Therefore the jet opening 
angle becomes narrow, which is determined by the in- 
verse of the Lorentz factor at the base of the free expan- 
sion, as explicitly shown by the numerical simulations. 
The opening angles are roughly constant over time with 
a factor ~ 2 fluctuations in Fig. [8[ 

We have also analytically analyzed the jet-breakout 
acceleration. The jet-breakout acceleration occurs with 
and without the collimating oblique shock, and both are 
equally important. For the former case, we solve the 
structure of the collimating oblique shock in a decreas- 
ing pressure profl le a nd obtain the Lorentz factor after 
the shock in Eq. (33). The post-shock Lorentz factor is 
enhanced appreciaEIy for a pressure slope close to A ~ 2 
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in P ex ^. The latter case happens after the last coUi- 
mating oblique s hock , even for A > 2 near the periphery 
of the jet inEq. (pSl). 



We have also compared our numerical results with the 
analytical fo rmulae for a jet propag ation inside the star 
presented by Bromberg et al. (2011 ), and have confirmed 
good agreements. I'br later use, we have calibrated the 
model parameter with the numerical results. We can now 
calculate the jet evolution relatively precisely with ease, 
such as the jet head position, the jet cylindrical radius, 
and the converging position of the coUimation shock, for 
a wide range of initial conditions. 

We have paid special attention to the numerical dif- 
fusion of the baryon loading into the jet through the 
discontinuity, which can entirely change the jet propa- 
gation. We have also taken the initial cylindrical radius 
of the jet sufhciently smaller than the to-be radius after 
the first coUimation shock because a large initial radius 
slows down the jet propagation. For these purposes, we 
have performed one of the highest resolution calculations 
so far. 

The post-breakout jet shows a hollow-cone angular 
structure. The edge is relatively sharp with an expo- 
nential drop. The bright rim is produced by the shock 
between the expanding jet and the high pressure cocoon 
before the free expansion. 

To understand the jet evolution, it is important to no- 
tice that the jet makes an off-center expansion as a result 
of the cocoon confinement. If the expansion origin is the 
stellar center, the fireball would feel a steep pressure pro- 
file of the cocoon pressure and the acceleration would be 
slow. In addition, the opening angle should be measured 
from the off-center origin for precise analyses. 

We have also applied our results to the observations. 



First we infer the initial Lorentz factor Fq of the jet in- 
jected at the central engine by using the observed open- 
ing angles in Fig. 



around 



2-3. 
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The distribution of Fq peaks at 
Second our result suggests the exis- 
tence of a maximum opening an gle for a high-entropy 
jet, 6'j^max ~ 1/5 ^ 12°, in Eq. (38). However several 
bursts violate this maximum value. This may imply a 
two-component jet with a baryon-rich slow sheath sur- 
rounding a central jet. 

We have also derived Amati and Yonetoku spectral re- 
lations by applying our results to the photosphere of the 
jet that is expanding freely after the jet breakout. We 
explain both the relations with only one assumption of 
a constant rate of the mass outflow M = Lj /Fq ^ const 



from the central engine in Eq. ( 42 1 . The fireball temper 



ature becomes different from the value at the central en- 
gine after the jet propagation through the star. Thus the 
confinement by the cocoon and the off-center expansion 
of the jet may be the missing pieces for the photosphere 
model so far. 

In the future, it is interesting to study the long-term 
evolution of the jet studied in this paper and the evo- 
lution of the low-luminosity jet and the two-component 
jet. It is also important to investigate the jet propaga- 
tion in a huge proge nitor such as the Populatio n III GRB 
( |Suwa fc Ioka||2011; Nagakura 
the ultra-long UKB ( |Levan et ai.||2013| ) 

We would like to thank A. Heger for his kindness 
to allow us to use his progenitor models. This work 
is supported by Grant-in-Aid for Scientific Research 
from the Ministry of Education, Culture, Sports, Sci- 
ence and Technology (MEXT) of Japan (20105005 (AM), 
24103006, 24000004, 22244030, 21684014 (KI)). 



Suwa & loka j2012|) and 



REFERENCES 



Abdo, A. A., et al., 2009, Science, 323, 1688 
Ackermann, M., et al., 2010, Astrophys. J., 716, 1178 
Aloy, M. A., Miiller, E., Ibanez, J. M., Marti, J. M., & 

MacFadyen, A. 2000, ApJ, 531, L119 
Amati, L., et al. 2002, A&A, 390, 81 
Begelman, M. C, & Cioffi, D. F. 1989, ApJ, 345, L21 
Bromberg, O., & Levinson, A. 2009, ApJ, 699, 1274 
Bromberg, O., Nakar, E., Piran, T., & Sari, R. 2011, ApJ, 740, 

100 
Campana, S., et al. 2006, Nature, 442, 1008 
Donat, R., & Marquina, A. 1996, Journal of Computational 

Physics, 125, 42 
Fong, W., Berger, E., Margutti, R., et al. 2012, ApJ, 756, 189 
Galama, T. J., et al. 1998, Nature, 395, 670 
Horiuchi, S., Murase, K., loka, K., & Meszaros, P. 2012, ApJ, 

753, 69 
loka, K., Murase, K., Toma, K., Nagataki, S., & Nakamura, T. 

2007, ApJ, 670, L77 
loka, K. 2010, Prog. Theor. Phys., 124, 667 
loka, K, Ohira, Y., Kawanaka, N., & Mizuta, A. 2011, 

Prog. Theo. Phys., 126, 555 
Iwamoto, K., et al. 1998, Nature, 395, 672 

Komissarov, S. S., & Falle, S. A. E. G. 1997, MNRAS, 288, 833 
Komissarov, S. S., & Falle, S. A. E. G. 1998, MNRAS, 297, 1087 
Lazzati, D., Morsony, B. J., & Begelman, M. C. 2009, ApJ, 700, 

L47 
Lazzati, D., Morsony, B. J., Margutti, R., & Begelman, M. C. 

2013, ApJ, 765, 103 

Levan, A. J., et al. 2013,"arXiv:1302.2352' 

MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262 

Marti, J. M. A., Mueller, E., Font, J. A., Ibanez, J. M. A., & 

Marquina, A. 1997, ApJ, 479, 151 



Matzner, C. D. 2003, MNRAS, 345, 575 

MazzaU, P. A., et al. 2006, Nature, 442, 1018 

Meszaros, P., & Waxman, E. 2001, Physical Review Letters, 87, 

171102 
Metzger, B. D., Giannios, D., & Horiuchi, S. 2011, MNRAS, 415, 

2495 
Mizuta, A., Yamada, S., & Takabe, H., 2004, ApJ, 606, 804 
Mizuta, A., Yamasaki, T., Nagataki, S., & Mineshige, S. 2006, 

ApJ, 651, 960 
Mizuta, A., & Aloy M. A. 2009, ApJ, 699, 1261 
Mizuta, A., Kino, M., & Nagakura, H. 2010, ApJ, 709, L83 
Mizuta, A., Nagataki, S., & Aoi, J. 2011, ApJ, 732, 26 
Morsony, B. J., Lazzati, D., & Begelman, M. C. 2007, ApJ, 665, 

569 
Murase, K., loka, K., Nagataki, S., & Nakamura, T. 2008, 

Phys. Rev. D, 78, 023005 
Nagakura, H., Ito, H., Kiuchi, K., & Yamada, S. 2011, ApJ, 731, 

80 
Nagakura, H., Suwa, Y., & loka, K. 2012, ApJ, 754, 85 
Pe'er, A., Ryde, F., Wijers, R. A. M. J., Meszaros, P., & 

Rees, M. J. 2007, ApJ, 664, LI 
Racusin, J. L., Liang, E. W., Burrows, D. N., et al. 2009, ApJ, 

698, 43 
Racusin, J. L., Dates, S. R., Schady, P., et al. 2011, ApJ, 738, 138 
Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 519, L17 
Shu, C.-W., &; Osher, S. 1988, Journal of Computational Physics, 

77, 439 
Stanek, K. Z., et al. 2003, ApJ, 591, L17 
Suzuki, A., & Shigeyama, T. 2013, ApJ, 764, L12 
Suwa, Y., & loka, K. 2011, ApJ, 726, 107 
Toma, K., loka, K., Sakamoto, T., & Nakamura, T. 2007, ApJ, 

659, 1420 



18 Mizuta & loka 

van Leer, B. 1977, Journal of Computational Physics, 23, 276 Yonetoku, D., Murakami, T., Nakamura, T., Yamazaki, R., Inoue, 

Wang, X.-Y., Razzaque, S., & Meszaros, P. 2008, ApJ, 677, 432 A. K., & loka, K. 2004, ApJ, 609, 935 

Woosley, S. E. 1993, ApJ, 405, 273 Zhang, W., Woosley, S. E., & MacFadyen, A. I. 2003, ApJ, 586, 

Woosley, S. E., & Heger, A., 2006, ApJ, 637, 914 356 

Zhang, W., Woosley, S. E., & Heger, A. 2004, ApJ, 608, 365 



